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Abstract. Weak gravitational lensing by the large scale structure can be used to probe the dark matter distribution 
in the Universe directly and thus to probe cosmological models. The recent detection of cosmic shear by several 
groups has demonstrated the feasibility of this new mode of observational cosmology. In the currently most 
extensive analysis of cosmic shear, it was found that the shear field contains unexpected modes, so-called B- 
modes, which are thought to be unaccountable for by lensing. B-modes can in principle be generated by an 
intrinsic alignment of galaxies from which the shear is measured, or may signify some remaining systematics in 
the data reduction and analysis. In this paper we show that B-modes in fact are produced by lensing itself. The effect 
comes about through the clustering of source galaxies, which in particular implies an angular separation-dependent 
clustering in redshift. After presenting the theory of the decomposition of a general shear field into E- and B- 
modes, we calculate their respective power spectra and correlation functions for a clustered source distribution. 
Numerical and analytical estimates of the relative strength of these two modes show that the resulting B-mode 
is very small on angular scales larger than a few arcminutes, but its relative contribution rises quickly towards 
smaller angular scales, with comparable power in both modes at a few arcseconds. The relevance of this effect 
with regard to the current cosmic shear surveys is discussed; it can not account for the apparent detection of a 
B-mode contribution on large angular scales in the cosmic shear analysis of van Waerbeke et al. (2002). 
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1. Introduction 

Gravitational lensing by the large-scale structure (LSS) 
leads to the distortion of the images of distant galaxies, 
owing to the tidal gravitational field of the matter inhomo¬ 
geneities. Following very early work on the study of light 
propagation in an inhomogeneous universe (e.g., Gunn 
1967; Kantowski 1969), Blandford et al. (1991), Miralda- 
Escude (1991) and Kaiser (1992) have pointed out that 
the observation of this ‘cosmic shear’ effect immediately 
yields information about the statistical properties of the 
LSS and, thus, on cosmology. Non-linear evolution of the 
matter spectrum was taken into account in later analytical 
(e.g., Jain & Seljak 1997; Bernardeau et al. 1997, Kaiser 
1998; Schneider et al. 1998, hereafter SvWJK) and numer¬ 
ical (e.g., van Waerbeke et al. 1999; Jain et al. 2000, White 
& Hu 2000) studies; see Mellier (1999) and Bartelmann & 
Schneider (2001; hereafter BSOl) for recent reviews. 
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It was only in 2000 when four teams nearly simulta¬ 
neously and independently announced the first detections 
of cosmic shear from wide-held imaging data (Bacon et 
al. 2000; Kaiser et al. 2000; van Waerbeke et al. 2000; 
Wittman et al. 2000). The detections reported in these 
papers (and in Maoli et al. 2001, using the VLT, and 
Rhodes et al. 2001, using HST images obtained with the 
WFPC2 camera) concerned various two-point statistics, 
like the shear dispersion in an aperture, or the shear 
correlation function. In van Waerbeke et al. (2001), the 
aforementioned statistics, as well as the aperture mass 
statistics (SvWJK), were inferred from the effective 6.5 
square degrees of high-quality imaging data. Very recently, 
Hammerle et al. (2002) reported on a cosmic shear de¬ 
tection using HST parallel images taken with the STIS 
instrument on an effective angular scale of ~ 30". 

The shear held, originating from the inhomoge¬ 
neous matter distribution, is a two-dimensional quantity, 
whereas the projected density held of the matter is a scalar 
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field. The relation between the shear 7 ( 0 ) = 7i(0)+i72(0) 
and the projected matter density k{6) is 


field K is related to the first spatial derivatives of the shear 
components in the following way (Kaiser 1995): 


7(6>) = - / d^e' V{e 

with the kernel 


e')n{9') , 


v{e) 


el-el- 216162 


( 1 ) 

( 2 ) 


here, k is the dimensionless surface mass density, i.e., the 
physical surface mass density divided by the ‘critical’ sur¬ 
face mass density, as usual in gravitational leasing; we fol¬ 
low the notation of BSOl in this paper. Since the two shear 
components originate from a single scalar field, they are 
related to each other; in particular, their partial deriva¬ 
tives should satisfy compatability relations, as we shall 
discuss in Sect. 2 below. In analogy with the polarization 
of the CMB, a shear field satisfying these compatability 
relations is called an E-mode shear field. 

Pen et al. (2002) pointed out that the cosmic shear 
data of van Waerbeke et al. (2001) contains not only an 
E-mode, but also a statistically significant B-mode con¬ 
tribution in addition. Such B-modes can be generated by 
effects unrelated to gravitational leasing, such as intrinsic 
alignment of galaxies (e.g.. Heavens et al. 2000, Crittenden 
et al. 2001a; Croft & Metzler 2000; Catelan et al. 2000) or 
remaining systematics in the data reduction and analysis. 

In this paper we show that a B-mode contribution to 
the cosmic shear is obtained by lensing itself. A B-mode 
is generated owing to the clustering properties of the faint 
galaxies from which the shear is measured. This spatial 
clustering implies an angular separation-dependent clus¬ 
tering in redshift, which is the origin not only of the B- 
mode of the shear, but also of an additional E-mode con¬ 
tribution. 

The paper is organized as follows: in Sect. 2 we provide 
a tutorial description of the E/B-mode decomposition of a 
shear field. Most of the results there were derived before in 
Crittenden et al. (2001b, hereafter COl), but we formulate 
them in standard lensing notation, which will be needed 
for the later investigation. The calculation of two-point 
cosmic shear statistics in the presence of source clustering 
is presented in Sect. 3 where it is shown that this clustering 
produces a B-mode. Numerical and analytical estimates of 
the amplitude of this B-mode are provided in Sect. 4 and 
discussed in Sect. 5. 


2. E/B-mode decomposition of a shear field 

In this section we provide the basic relations for the de¬ 
composition of the shear field into E- and B-modes. Most 
of these relations have been obtained in COl; we shall write 
them here in standard lensing notation. 

2.1. Motivation 

If the shear field is obtained from a projected surface mass 
density k as in Eq. (|^), then the gradient of the density 


Vk = 


7 1.1 + 72,2 

72.1 — 7l,2 


= u 


(3) 


The vector field u can be obtained from observations, e.g. 
in weak lensing cluster mass reconstructions, by obtain¬ 
ing a smoothed version of the shear field and then differ¬ 
entiating this numerically. Owing to noise, the resulting 
(‘observed’) field u will in general not be a gradient field. 
The non-gradient part of u is then a readily identifiable 
noise component and can be filtered out in the mass re¬ 
construction. Seitz & Schneider (1996) provided a scheme 
for this noise filtering (see also Seitz & Schneider 2001 for 
a simpler though equivalent method), which was shown by 
Lombardi & Bertin (1998) to be an optimal reconstruction 
method. 

If the shear field cannot be ascribed to a single geomet¬ 
rically thin gravitational lens, the non-gradient part of u 
is not necessarily due to noise. For example, if the galaxies 
have intrinsic alignments, this may induce a curl-part of 
u. To project out the gradient and curl part of u, we take 
a further derivative of u, and define 

= V • u ; = V X u = U 2 ,i — Ml, 2 ■ (4) 

Through these relations, k® and kP are not uniquely 
defined on a finite data field; as discussed in Seitz & 
Schneider (1996), a further condition is needed to specify 
the two modes uniquely. However, we shall not be con¬ 
cerned here with finite-field effects. 

An alternative way to define k® and is through the 
Kaiser & Squires (1993) mass-reconstruction relation 


K^{e)+iK.^{9) = - [ d'^e'V*{9-9')-f{9') , (5) 


which formally requires data on an infinite field; here, T>* 
denotes the complex-conjugate of the complex kernel (||) . 
If 7 is of the form ( 0 ) with a real field k, then the result 
from (||) will be real, k® = k, = 0. In applications 
of the KS-formula i) to observational data, where the 
recovered shear field necessarily is noisy, one usually takes 
the real part of the integral to obtain the projected mass 
density field. For a general shear field, the result from 
will be complex, with the real part yielding the E-mode, 
and the imaginary part corresponding to the B-mode. 

To simplify notation and calculations, it is convenient 
to express two-component quantities in terms of complex 
numbers. We define the E- and B-mode potentials and 
V’® by 

VV®’® = 2«:^’® , (6) 


and combine the two modes into the complex fields 

K = + Ik^ , ip = + iV'^ • (7) 


The complex shear 7 = 71 -I- 172 is obtained from the 
potential ■0 by 7 = Dip, where the differential operator 
D = (9ii - 922)72 -I- i9i2; hence, 


7 = 


- P%) - 


+ i 


^:^i2+^ 
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Inserting this into i) yields 


u = 


+ '‘\1 ) 


Indeed, the shear field can be decomposed into E/B- 
modes, 7 = 7 ^ + i 7 ®, with 


p 'D 
7 = — *TZe 

TT 


V* 

— * 7 
TT 


R 2 ? 

7 ^ = — * Xm 

TT 


V* 

— * 7 
TT 


where the operator T> is defined in Eqs. ( 0 ), (|), and 
denotes complex conjugation. Thus, the two components 
can be obtained from the shear field by filtering, except 
for an additive constant. 


2.2. Shear correlation functions and power spectra 

The discussion above dealt with the shear field itself. In 
the application to cosmic shear, one usually does not in¬ 
vestigate the shear of a K-field itself, but its statistical 
properties. In this paper we shall concentrate solely on 
two-point statistical measures of the cosmic shear, and 
their decomposition into E- and B-modes. 

Owing to statistical homogeneity and isotropy of the 
Universe, are homogeneous and isotropic random 

fields. Hence, in terms of their Fourier transforms 

= J , ( 8 ) 

one defines the two power spectra Pp, Pb, and the cross 
power spectrum PpB by 

= i2Tr)^ dn{e - £') PEii), 

{k^{£)k^*{£')) = {2Tr f Su{i - t) Pb w , (9) 

{k^{£)k^*{£')) = {2^f5E{t-t)PEn{f), 

where fo denotes Dirac’s delta distribution. In terms of 
the complex field k, we then have 

{k{£)k*{£')) = {2nf5E{f-f')[PE{f) + Piim , 
{k{£)k{£')) = {2nf5E{e. + f') ( 10 ) 

X [PE(^)-PB(^) + 2iPEB(^)] . 

The Fourier transform 7 (f) of the shear is related to k{£) 
through 




( 11 ) 


where (3 is the polar angle of £. The correlators of the 
shear then become 


{l{£)T{f')) = {2nf5E{i-l')[PE{f) + PB{(.)\ , 
imW)) = (27r)2<5D(^ + /)e4'^ (12) 

X [Pe(^) — Pb(^) + 2iPEB(^)] • 


Next we define the correlation functions of the shear. This 
is done by considering pairs of positions and 0-1-19, and 
defining the tangential and cross-component of the shear 
at position i9 for this pair as 

7 t =-^ 0 ( 70 “^“^) , 7x =-2:m(7e“^“^) , (13) 

respectively, where (p is the polar angle of the separation 
vector 6. Then, the shear correlation functions are defined 
as 


^+{^) = (7t7t) + (7x7x) {0) , 

C-(^') = (7t7t) - (7x7x) {8) , (14) 

fx{8) = (7t7x) {9) ■ 


The shear correlation functions are most easily calculated 
by choosing 6 = {6, 0 ), in which case 7 t = — 71 , 7 x = —I 2 , 
and expressing the shear in terms of its Fourier modes. 


(7(O)7*(0)) = 

d^£ f d2f Q 




(27r)2 J (27r)2 
d££ 

hm[PE{£) + PB{e)] , 


Jo 27 r 

(7(O)7(0)) = e-(0) + 2iex(0) 

f d^£ r d^e 


i£'e 


J (2WJ (2W 7 (<) 7 («')> ( 16 ) 

d££ 

Jo ^ ~ ^ 2iPEB(^)] . 


Making use of the orthogonality of Bessel functions, 
[ d9 e3^{s0)3^{t9) = ^ 


(17) 


we can invert the relations (|T^) and ( |l^ and express the 
power spectra in terms of the correlation functions, 

pOO 

PE{i) = 7 t deo , 

Jo 

poo 

PbW = tt / [^^{9)Jom - . (18) 

pOO 

Peb(£) = 27r / d9 9f^i9)34m, 

Jo 


Hence, we have now expressed the various power spectra 
in terms of the directly observable correlation functions 
f. One notes that the correlation functions and de¬ 
pend on both, the E- and B-mode power spectra, whereas 
the cross-correlation ^x depends on the cross-power Peb 
only. It is obvious that the cross-power and its correspond¬ 
ing correlation function do not ‘mix in’ with the E- and 
B-mode; in addition, the cross-power vanishes if the shear 
field is statistically invariant under parity transformations, 
which leave 7 t unchanged, but transform 7 x ^ — 7 x • One 
can therefore assume that in realistic cases, ^x = 0 = Peb- 
However, since cosmic shear is measured from finite data 
fields, cosmic variance may lead to a non-zero measure¬ 
ment of the cross-power; in fact, the measurement of the 












4 


Peter Schneider et al.: B-modes in cosmic shear from source redshift clustering 


cross-power may serve as a lower limit on error bars of the 
other power spectra. 

For a determination of the power spectra, the expres¬ 
sions require a measurement of the correlation func¬ 
tions over an infinite range in angle; whereas the correla¬ 
tion functions decrease with 0 and become very small for 
large 9, so that in effect the integrals can be replaced by 
ones over a finite range of integration, one might want to 
obtain more local decompositions into E- and B-modes. 


2.3. E/B-mode correlation functions 

We define the four correlation functions 
f 

fE.B+iO) = J PE,Bif) doidf) , 
f dll 

f,E.B-{0) = J PE,Bil) diiOl) , (19) 

which are defined such that in the absence of B-modes, 
^E± = C±; these four correlation functions have also been 
defined in COl, although only the ‘-I-’ ones were investi¬ 
gated in more detail there. Inserting ( 0 ) into the forego¬ 
ing definitions, one obtains 

-I pOO pOC 

fE+{9) = - dll ddd3o{9l) 

^ Jo Jo 

X [5+(r?)Jo(d£)+e-(t?)J4(^?^)] . (20) 


The ^-integration can be carried out; consider the function 

pOO 

G{'d,9)= / dt 13o{td) 34{t6) . (21) 

JQ 

Making use of the recurrence relations for Bessel functions, 
one can express J 4 as 

24 8 

J 4 (a:) = -^32{x) -Ji(a;) -I- Jo(a;) . 

X 


By using Eq. (11.4.41) of Abramowitz & Stegun (1965), 
together with 0 ), one can perform the integration in 
( 0 ) term by term to obtain 

/4 12i9^\ 1 

^ j ’ ( 22 ) 

where H(a;) is the Heaviside step function. We also note 
the interesting property, 


pOC 

/ dddG{d,9)G{0,(p) = 

Jo 


Su{0 - if) 




= / dddG{9,d)G{(p,'0), (23) 

Jo 

which is readily shown using (0)- Thus, becomes 

1 


fE+if^) = 


M0) 


dd0f-id)G{9,d) 


poo in / n 2 

M0) + U0) + l ^Ud)(d-i2- 


.(24) 


order to obtain ^e-i- one would need to know for arbi¬ 
trarily large separations. However, 


fE-iO) = 


1 


UO) 


d0df+i'd)Gi'd,9) 


U0) + ^+io) 


ddd 


?+(^^) 4-12 




(25) 


depends on the observable correlation functions over 
a finite range only and thus can be measured from hnite 
data sets. Analogously, we hnd for the B-mode correlation 
functions 


eB+(0) = i ?+(0)- pd0df_{i})G{e,0) 

poo 

-i-{e)+ / dddf+{d)G{d,e) 

Jo 

f'^ ddd ( dP 

C+(0)-e-(0)+ I -^M0)(d-12- 


fB-id) = 


(26) 


We note that ^e-i- + Cb+ = Cs-i Ce- — Cb- = f-- In order 
to calculate the E- and B-mode correlation functions, one 
needs to know either the observable correlation function 
f- to arbitrarily large, or to arbitrarily small sepa¬ 
rations. This is of course impossible, owing to the hnite 
size of data helds on the one hand, and the impossibility 
to measure shapes of very close pairs of galaxies. In ei¬ 
ther case, the lack of measurements for large (or small) 
separations can be summarized in two constants: suppose 
that can be measured down to separations of 6 *i„in; 
then, the integral in the ’ modes over can be split 
into one from 0 to 0min, and one from 0niin to 9. The for¬ 
mer one has the 0 -dependence a/9‘^ — 9/9^^, where a, b are 
two constants, depending on for 0 < 0min. The decline 
of this contribution, with leading order 0 “^, shows that 
it has a small inhuence on the determination of the ’ 
modes for 0 0min; in addition, ‘reasonable guesses’ for 
a and b may be obtained by extrapolating the measured 
towards small angles. The same reasoning shows the 
analogous situation for the ‘-I-’ modes. 


2.4. Aperture measures 

One very convenient way to separate E- and B-modes is 
provided by the aperture mass: Defining the tangential 
and cross component of the shear relative to the center of 
a circular aperture of angular radius 0 , and defining 

Map(0) = J d^dQi\d\)j,{d) , 

Mj.( 0 ) = J d20Q(|,9|)7x(i9) , (27) 


We have obtained a combination of shear correlation func- it was shown in COl that E-modes do not contribute to 
tions which depends only on the E-modes; however, in M±, and B-modes do not contribute to Map; in fact, this 
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can be easily seen directly by inserting the Fourier trans¬ 
form of the shear © into (p^). Here, Q{'&) is an axially- 
symmetric weight function which can be chosen arbitrar¬ 
ily. The integration range in the foregoing equations ex¬ 
tends over the support of the weight function Q. The aper¬ 
ture mass was introduced by Schneider (1996) in an at¬ 
tempt to detect mass peaks from shear fields, and later 
used by SvWJK as a two- and three-point statistics for 
cosmic shear. Map can also be written as a filtered version 
of the surface mass density (with a different and compen¬ 
sated weight function which is related to Q), whereas M± 
has no direct physical interpretation; in the absence of B- 
modes, M± should vanish, and any non-vanishing signal 
is usually interpreted as being due to noise or remaining 
systematics, and thus as a convenient error estimate for 
Map. In SvWJK, a family of convenient weight functions 
Q was considered, the simplest of which is 

where 9 is the radius of the aperture. This form of the 
weight function shall be assumed in the following. 

Using the complex number 

M(d) = Map(d)+iM^(d) 

= -J d2dQ(|r?|)7(i9)e-^‘‘", 


where (p is the polar angle of i?, one finds that 


(Mfp) + (Mi) = (MM*) 

= ^ [ ^ [Pe(^) + -Pb(^)] W{ei) , (28) 

Jo 

with 


W{v) := 


576J1(?7) 

-n4 


(29) 


which was derived by using the Fourier transform ©) of 
the shear, and the final steps are as in SvWJK. Similarly, 
one obtains 


(Mfp) - (Mi) + 2i (MapMi) = (MM) 

= ^ / de£ [Pe(^) - Ub(^) + 2iPEB(^)] wm . (30) 

Jtt Jq 

Combining the two previous equations, one thus gets 


1 

(Ma"p)(d) = —y de£PEie)wm, 

(Mi)(d) = — I de£PB{£)Wi0£) , (31) 


so that these two-point statistics clearly separate E- and 
B-modes. In addition, they provide a highly localized mea¬ 
sure of the corresponding power spectra, since the filter 
function W{ij) involved is quite narrow (see SvWJK); in 
fact, Bartelmann & Schneider (1999) have shown that re¬ 
placing W{rj) in ( ^ ) by ASuir] — r]o), with r]o « 4.25, pro¬ 
vides a fairly accurate approximation. Furthermore, (^o|) 


can be used to check whether the shear data contain a 
contribution from the cross power Peb- 

The aperture measures can be obtained directly from 
the observational data by laying down a grid of points, 
at each of which Map and M_l are calculated from (^. 
However, obtaining the dispersion with this strategy turns 
out to be difficult in practice, since data fields usually con¬ 
tain holes and gaps, e.g. because of masking (for bright 
stars), bad columns etc. It is therefore interesting to cal¬ 
culate these dispersions directly in terms of the correlation 
functions, which can be done by inserting (^^ into 




(Ml) (9) = ^ 


dd-d 

02 

ddi? 






+ ^-{d)T. 
-C-WT- 



(32) 


where we have defined the functions 

/»00 1 I 

T+{x) = 576 J ■^Jo(M) [J 4 (t)]^ , 

/•OO 1 . 

r_(x) = 576 J [J4(t)]^ ■ 


(33) 


The integration range in (32) formally extends from zero 
to infinity, but as we shall see shortly, the functions T±{x) 
vanish for x > 2, so the integration range is 0 < < 29: 

For P_(x), an analytic expression can be obtained, using 
Eq. (6.578.9) of Gradshteyn & Ryzhik (1980), 


r_(x) = —x^ 

^ ^ 357r 


1 - 


7/2 


H(2 - x) , 


(34) 


so that P-(x) vanishes for x > 2. Eurthermore, the two 
functions T+ and P-(x) are related: using (|r^ , one finds 
that 

nOC poo 

/ dx X P+(x) Jo(tx) = IU(t) = / dx xP_(x) J 4 (tx) 

Jo Jo 

so that 

poo poo 

r_(x) = / dttw{t)34{xt) = dy yT+{y)G{y,x) , 

Jo Jo 

poo poo 

T+{x) = / dttW(t)3o{xt) = dy yT_{y)G{x,y) ; 

Jo Jo 


using the latter expression, together with 
one obtains 


and 


P+(x) = 


6(2-15x2) 


1 -arcsin 

TT 


(I) 


c\/4 — a 
IOOtt 


(120 -f 2320x2 - 754x‘‘ -h 132x® - 9x®) (35) 


for X < 2, and T+(x) vanishes for x > 2. Hence, the inte¬ 
grals in extend only over 0 < < 20, so that (Mfp) 

and (m2 ) can be obtained directly in terms of the ob¬ 
servable correlation function over a finite interval. The 
two functions T± are plotted in Fig. 1. 
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Fig. 1. The four functions defined in text 


2.5. Shear dispersion 

Another cosmic shear statistics often employed is the shear 
dispersion in a circle of angular radius 6. It is related to 
the power spectra by 

(0) = ^ [ d« (Pe + Pb)W Wrnm , (36) 


-|2 


where 


^th(^) = 


27T 


43lir,) 


(37) 


is the top-hat filter function. In contrast to the aperture 
measures of the previous subsection, the shear dispersion 
( ^ ) contains both modes; furthermore, the filter function 
W'TYii'n) is much broader than W{rf) in (^), as demon¬ 
strated in SvWJK. It thus provides a much less localized 
measure of the power spectra than the aperture measures. 
On the other hand, this larger filter width implies that 
the signal of the shear dispersion is larger than that of 
the aperture measures, which explains why the first cos¬ 
mic shear detections (van Waerbeke et al. 2000; Bacon et 
al. 2001, Kaiser et al. 2000) were obtained in terms of the 
shear dispersion. 

As before, the shear dispersion can be obtained by cal¬ 
culating the mean shear in circles which are laid down on 
a grid of points, with the drawback of being affected by 
gaps in the data field. Alternatively, the shear dispersion 
can be obtained directly from the correlation function. 




<"' = / 


6»2 ( 0 


where (van Waerbeke 2000) 

S+{x) = — 4arccos f— x \/4 — : 
TT L V 2 / 


(38) 


(39) 


for X <2, and zero otherwise. Hence, the integral in (38) 
extends only over the finite interval 0 < r? < 20, which 
makes this a convenient way to calculate the shear disper¬ 
sion. 

One can also define the shear dispersions of the E- and 
B-mode, according to 






(40) 


but they cannot be individually obtained from measuring 
the shear directly. Nevertheless, both of these dispersions 
can be obtained in terms of the correlation functions. 


ItI" 


/■( 

{ 0 )= ■ 
E,B J 


'W 




±^-{d)s. 


(41) 


which can be derived in close analogy to the derivation 
of (1^), and the function S- is related to S+ in the same 
way as the corresponding T-functions, 


pOO 

S-{x) = / dy yS+{y)G{y,x) 
Jo 


x^jA — a;2(6 — x"^) — 8(3 — a;^) arcsin(a:/2) 


(42) 


for a; < 2, and ^-(a;) = 4(a;^ — 3)/a:^ for x > 2. Hence, 
the integrals in (^ ) do not cut off at finite separation, 
which was to be expected, since a constant shear cannot 
be uniquely assigned to an E- or B-mode, but contributes 



3. B-mode from source clustering 

In the previous section we have presented the decomposi¬ 
tion of a general shear field into E/B-modes. It is usu¬ 
ally assumed that lensing alone yields a pure E-mode 
shear field, so that the detection of a B-mode in the van 
Waerbeke et al. (2001) data (see also Pen et al. 2002) was 
surprising and interpreted as being due to systematic er¬ 
rors or a signature of intrinsic alignment of sources. Here 
we show that lensing indeed does generate a B-mode com¬ 
ponent of the shear if the source galaxies from which the 
shear is measured are clustered. 


3.1. Correlation functions and power spectra 


Define the ‘equivalent’ surface mass density for a fixed 
source redshift, or comoving distance w, 

pW 

k{6,w) = / dw' F{w',w) S[f{w')6,w'] . (43) 

Jo 


where 
F{w', w) 


3ggHo f{w')f(w - w') ^ 
2 c2 a{w')f{w) ’ 


(44) 


here, Hq and Hg denote the Hubble constant and the den¬ 
sity parameter, w is the comoving distance, f{w) is the 





















Peter Schneider et al.: B-modes in cosmic shear from source redshift clustering 


7 


comoving angular-diameter distance to comoving distance 
w, S is the density contrast, and a{w) = (1 -|- z)~^ is the 
cosmic scale factor, defined such that a = 1 today, again 
using the notation of BSOl. Accordingly, we define the 
shear components 

7a(0,w) = (0,ui) , (45) 

where the operator V is defined in (|^ and (^). Then, the 
shear correlation function for two sources at positions 0 i 
and distances Wi becomes 

( T) T)a \ 

“2 ) * J Piw[,wi) 

pW 2 

X dW 2 F{w 2 ,W 2 ) {S[f{w[)ei,w[] 5 [/(u;2 )02,W2]) , 

Jo 


Wi seen from the distance w < Wi and that seen from an 
observer at w = 0. 

When measuring cosmic shear from source elliptici- 
ties, the source galaxies have a broad distribution in red- 
shift, unless information on the redshifts are available and 
taken into account. Hence, to calculate the observable 
shear correlation functions, the foregoing expressions need 
to be averaged over the source redshift distribution. Let 
p{wi,W 2 ', 0 ) be the probability density for comoving dis¬ 
tances of two sources separated by an angle 0 on the sky; 
then we have for the observable correlation functions 

C±(6') = / dwi du ;2 p(i«i, ((7t7t)±(7x7x)) (48) 

Jo Jo 

and wh is the comoving distance to the horizon. By chang¬ 
ing the order of integration according to 


where the first V operates on 0i, and the second V on 
62 - This equation is of the form (BS 2.78), and thus we 
obtain, using (BS 2.82), 

VaVn /■“ 1.2 

ilalp) = ^ * / dw F{w,Wl) F{w,W 2 ) 

^ Jo 

/ d'^k 

P 5 (|k|,w)exp [i/(w)k • {Oi - 62 )] ,(46) 

where we have temporarily dropped the arguments of the 
shear correlator, wi ^2 = min(r(;i,^ 2 ), and Ps is the power 
spectrum of the density fluctuations which develops as a 
function of cosmic time (or as a function of comoving dis¬ 
tance w). The upper limit of the integral expresses the fact 
that a correlated shear can only be generated by matter 
which is at smaller distance than both sources. 

The operators V only act on the final term of ( ^^ 
which can be evaluated using the Fourier transform of V, 
as in (ITlI), 



we obtain 


U{ 0 ) 


X 


dw r ^ / i 

40"^ Jo J (27r) \f{w) 

Jo,4(^^) (^(^^7 a^i) ^(^a,'*^ 2 )) 7 



(49) 


where the angular brackets denote the averaging of the 
angular-diameter distance ratios over the source distance 
distribution. We shall write the source redshift distribu¬ 
tion as 

1 -I- UJ{0) 

where 

(^{0) = / dw p^{w) g{w;9) (51) 

^0 


so that 

dui /(wi — w) f{w 2 — w) 




4c4 Vo a‘^{w) /(mi) f{w 2 ) 

Evaluating the relevant combinations, one finds 

// X / XI./, X 9^0dm 

(( 7 t 7 t) + ( 7 x 7 x)) (^';'a;i,m2) = ' 


4c'* Jq a?‘{w) 

r d££ / £ \ 

xR{w,wi) R{w,W 2 ) / TTTT Ps (-rr^,w] do{£0), 


(27r) V/M’ 


((7t7t) - (7x7x)) {0;wi,W2) = 


dm 


4c'* 


\w) 


is the angular correlation function of the galaxies, and 
p(m) describes their redshift distribution. The second term 
in ( ^ accounts for source clustering. In making this 
ansatz, we have accounted for the fact that redshift clus¬ 
tering occurs only over a very small interval in redshift 
over which all the other redshift-dependent functions oc¬ 
curring in ( ^ can be considered constant. Note that (M) 
is normalized, 


ydmi J dw2 p{wi,W2-,9) = I , 


as required. Then, the average of the angular-diameter 
distance ratios becomes 


f d£i ft \ 

xi?(m, mi) i?(m, m2) J ^j J4(£6'), 

where 9 = \6i — 02| and R{w,Wi) = f(wi — w)/ f(wi) is 
the ratio of the angular diameter distances of a source at 


(i?(m, mi) i?(m, m 2 )) = 
where 

rwn 

W{w)= / dmi p(mi) 


W‘^{w) + V{w,9) 


1 -I- uj{9) 

/(mi - m) 

f{wi) 


(52) 


(53) 


pwn 

V{w,9) = / dmi p^(mi) 

J W 


/(mi - m) 
fiwi) 


giwi-9) . (54) 
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The correlation of sources thus yields an average of the 
angular-diameter distance ratios which is not simply the 
square of the mean distance ratio W, but contains in ad¬ 
dition a correlated part described by V and the normal¬ 
ization correction 1 -|- w. If the angular separation of the 
sources is large, the correlation in redshift is expected to 
be small; hence, for large separations one expects V and 
(jj to vanish. The degree of redshift correlation depends on 
the angular separation considered; the fact that the mean 
of the product of the angular diameter distance ratio 
depends on the separation 0 is the cause for a B-mode 
contribution to the shear correlation function! 

One can check that the correlated redshift probabil¬ 
ity distribution behaves as expected in some simple cases. 
For example, if w c/Hq is very much smaller than the 
characteristic source distance wg, one finds that 

{R{w,wi) R{w,W 2 )) ~ 1 ; 


in this case, the leasing strength of matter at distance 
w is basically independent of the exact source redshift, 
so that source redshift clustering is irrelevant for those 
lens redshifts. Another case of interest occurs when the 
selected sources come from a very narrow distance inter¬ 
val, of width Aw centered on wg; then, d) yields the 
relation uj{0) « giwg; 9)/Aw, and 


{R{w, wi) R{w, W 2 )) 


f{wg - w) 
f{wg) 


H(wo - w) . 


Hence, also in this case, (RR) does not depend on 9, and 
therefore no B-mode contribution occurs - as noted before, 
if all sources are at the same redshift, one obtains a pure 
E-mode shear field. 

We can now rewrite (^) in the form 

1 f d£ £ 

= TT^ J [P.ii) + Pcii;0)] , (55) 

where the 0 (4) corresponds to (C-)) 


P.{£) 


40“^ 



dw 

a^(w) 


W^(w)Ps 



(56) 


and 

(tr'”) 

The first term of ( |^ in the absence of source correlations 
(i.e., w = 0) is the one usually derived in cosmic shear 
considerations; is the power spectrum of the projected 
matter density, related to the three-dimensional power 
spectrum Ps by a Limber-type equation (e.g., Kaiser 
1992). The second term in (^ ) and the ‘normalization 
correction’ 1 -|- w comes about due to source correlations. 


we can derive the E- and B-mode power spectra, making 
use of (|l^), 

Pe(£) = P.(£) + Pce(£) ; Pb(£) = Pcb(£) ; (59) 

with 

PcE,bW = I 

X [Jg{£9)dgi£'9) ± J4{£9)34{£'9)] , (60) 


where the ‘-I-’ (‘—’)-sign corresponds to the E- (B-) mode. 
If the ratio containing the power spectra did not depend 
on 9, Pb would vanish identically, as it should. However, 
this term does depend on 9 due to source correlations; 
therefore, without using redshift information, the presence 
of B-modes in cosmic shear observations is unavoidable. 

Using the definitions of the E- and B-mode correlation 
function, we obtain 


^E+{(^) = 


-k 


d££ 

(27r) 

d££ 

(dir) 


P.{£) 3om 

P{£-9) [3oi£9) + 34i£9)] 


A -| 0^2 

dd£}Pi£;£))34m[^-—, 


^B+i(^) = 


d££ 


V{£-,9) [3g{£9)-34{£9)] 

4 1202 


(471-) 

X3 

d£)£)P{£-,d)34{£d) 


02 04 


^E-{0) = 


d££ 

(27r) 

d££ 


P.{£)d4{£0) 

V{£- 9) [3g{£9)+34m] 


J ( 47 r) 

^"d00P(£;0)Jo(£0) (^^- 


120^ 

~9P 


Cb-(6') = 


d££ 

(47r) 


P{£-,9) [3gi£0)-34m] 


^'d00P(^;0)Jo(e)(^^- 



1202 y 

V 02 

04 ) 


(61) 


(62) 


(63) 


(64) 


3.2. The E/B-mode decomposition 

Erom the correlation functions (|^), by writing 
P^{£) + Pc{£;9) 


1 -k U!{9) 


P,{£;9) - w{9)P,{£) 


= P.i£) + 


1 -k uj{9) 


4. Relative strength of the B-mode 

Whereas the presence of a B-mode, and an additional con¬ 
tribution to the E-mode due to source clustering must oc¬ 
cur, one needs to estimate the relative amplitude of this 
effect as compared to the ‘usual’ cosmic shear strength 
described by This estimate requires a model for the 
source clustering, i.e., a model for the function g{w,9). g 
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can be related to the three-dimensional correlation func¬ 
tion ^gg of galaxies, 

g{w-e) = 1“ d(Au;) ^g + p{w)e^) . (65) 

If we assume that the correlation function behaves like a 
power-law, i^g(r) = [r/ro(w)]~^, where rpw) is the co¬ 
moving correlation length, and 7 ~ 1.7, then 

g{w-,e)=Cro{w)(^^^^^ , (66) 

where C = v^r ([7 — l]/ 2 )/r( 7 / 2 ), and r(a:) is the 
Gamma function. This yields for the angular two-point 
correlation function 

ui{0) = [ dw p'^{w)rQ{w) ( ^ . (67) 

Jo \mw)J 

A useful parameterization of oj{9) is w = A(l')(6*/1')(^“^). 
Fixing 7 = 1.7, one obtains a relation between the corre¬ 
lation length ro(iy) and the redshift distribution p(w). 



1 


dw p^{w) ro{w) 


f fjw) 

\ro{w) 


- 0.7 


= 1.65x10 


-5 -4(10 
0.02 


where the fiducial value of A(l') was taken from 
McCracken et al. (2001). Note that in McCracken et 
al. ( 2001 ), essentially the same data set has been used as in 
the cosmic shear analysis of van Waerbeke et al. (2001); 
in particular, the depth of the data are the same. The 
above-quoted value for the angular clustering strength at 
1 ' corresponds to the faintest flux threshold considered in 
McCracken et ah, which is very similar to the flux limit 
employed in the cosmic shear analysis. It must be men¬ 
tioned, however, that the galaxies used in the cosmic shear 
analysis do not form a truly flux-limited sample, since ad¬ 
ditional cuts are used, e.g. a size cut. Hence, a precise esti¬ 
mate of the angular correlation function of those galaxies 
which were used for the cosmic shear analysis cannot be 
given. 

For this reason, we shall assume the power-law depen¬ 
dence of io{0) as given above; in addition, we will make 
the simplifying assumption that the comoving clustering 
length rQ{w) is independent of distance w] this assump¬ 
tion is not too critical, since the function p^(w) is rela¬ 
tively well peaked and therefore large ui-variations of the 
correlation length are not probed. Then, ( | 6 ^ determines 
this constant comoving correlation length rp. We obtain 
in this case 


Fig. 2. Dimensionless power spectra PP{£), as a func- 
’ (*^^iion of of wavenumber i. The solid curve corresponds to 
the power spectrum Ppptj that is the ‘standard’ power 
spectrum of the projected mass density. The dotted curve 
displays PPc{£, l')j and the two dashed curves correspond 
to the E- and B-mode power caused by the source cluster¬ 
ing. Here, a ACDM model was used, with shape parameter 
F = 0.21, normalization erg = 1, and the source redshift 
distribution is characterized by zq = 1, yielding (z) « 1.5. 
Other parameters for the model used here are mentioned 
in the text. 


g{w;e) = uj{e) 




1-7 


/ dw' [/(ru)]^” 


the 0 -integration can then be carried out first, making use 
of Eqs. (11.4.33, 34) of AS. 

In order to make further progress, we need to assume a 
redshift distribution for the sources from which the shear 
is measured. We employ the form (Brainerd et al. 1996) 


Pz(z) cx z exp 


zo 


(69) 


The power-law dependence of g on 0 implies that Pc{£;0) 
also behaves like Since the angular correlation func¬ 

tion is small compared to unity, even on scales of a few 
arcseconds, we shall neglect uj{6 ) in the denominator of 
the integrand in (60); this greatly simplifies the calcula¬ 
tion of the power spectra due to source clustering, since 


and shall consider /3 = 3/2 in the following, yielding a 
mean redshift of (z) « 1.5zo. 

In Fig. 2 we show an example of the various power 
spectra considered here; all power spectra are multiplied 
by P to obtain dimensionless quantities. For this figure, 
we employed a standard ACDM model with Ha = 0.7 
and Hq = 0.3 and normalization erg = 1. Sources are dis¬ 
tributed in redshift according to the foregoing prescrip¬ 
tion, with zq = 1. The amplitude of the angular correla¬ 
tion function of galaxies was chosen to be A(l') = 0 . 02 , 
and the slope of 7 = 1.7 for the three-dimensional cor¬ 
relation function was used; the corresponding correlation 
length in this case is tq « 4.7/i“^Mpc. To calculate the 
three-dimensional power spectrum and its redshift evolu- 
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tion, we used the Peacock & Dodds (1996) prescription 
for the non-linear evolution of Ps{k,w). 

The power spectrum of the projected mass density, 
Pk{^), is the same as that in SvWJK, except for a slightly 
different choice of the cosmological parameters. The spec¬ 
trum Pc{i', 10 is very much smaller than Pk(£), as expected 
from the smallness of the amplitude of the angular corre¬ 
lation function; in fact, the ratio T’c(^; lO/-fK(^) i® nearly 
constant at a value of approximately (1 -I- B)A{1'), with 
B « 1.2 for this choice of the parameters. 

The behavior of the power spectra which arise from 
source clustering, PcE and Pb = PcB, as a function of i is 
quite different. First, both of these spectra are very sim¬ 
ilar, which is due to the fact that the J 4 -term in (|6^ ) is 
much smaller than the Jo-term. Second, although both of 
these spectra are small on large angular scales, i.e. at small 
their relative value increases strongly for smaller angu¬ 
lar scales. Hence, as expected, the relative importance of 
source clustering increases for larger £. What is surpris¬ 
ing, though, is that these power spectra have the same 
amplitude as P^ at a value of £ ~ 6.7 x 10^, corresponding 
to an angular scale of 0 = 211 ji ^ 2", and the relative 
contribution of the B-mode amounts to about 2% at an 
angular scale of 1'. It should be noted here that cosmic 
shear has already been measured on scales below 1'; there¬ 
fore, source clustering gives rise to a B-mode component 
in cosmic shear which is observable. 

We shall now consider the behavior of Pb{£) for large 
values of i. The aforementioned properties of Pc{£', 0) can 
be summarized as 

Pc(£;0)«(l + P)H(1')P«(£)(^) \ 

Inserting this result into (^), neglecting the J 4 -terms and 
considering the limit ^ oo yields 

fPB{£)^aiBA{l') J diipyi) , 


where Oi = 2i-^r[{3 - 7)/2]/r[(7 - I)/2] w 0.335 for 
7 = 1.7. Hence, Pb oc at large £. To obtain an ap¬ 
proximate value for the integral in the preceding equation, 
we shall describe the power spectrum P^ by a simple func¬ 
tion, 


fpye) 




with Pm ~ 3.7 X 10 a ~ 0.7, [3 ~ 0.6, and £m ~ 7 x 10^. 
Using H(l') = 0.02, this then yields 


fPB{£) 


1.5 X 10""^ 



which is a reasonably good description of the result in 
Fig. 2 for large i. Furthermore, we can obtain the ratio 
Pb/Pk in the limit of large f » £m, which yields 

Pb(^) ^ / I Y " 

P,,(£) V4-9 X lOV ’ 



Fig. 3. For the same model as in Fig. 2, several correlation 
functions are plotted. The solid line shows C-i-(^); in fact, 
the correlation function ^e-i- cannot be distinguished from 
^_i_ on the scale of this figure; their fractional difference 
is less than 1%, even on the smallest scale shown. The 
two B-mode correlation functions are shown as well as 
and ^E-- Note that the difference between the latter two 
is larger than that of the corresponding ‘-I-’-correlation 
functions. 


and roughly predicts the correct crossing point between 
these two power spectra seen in Fig. 2. 

In Fig. 3 we have plotted several correlation functions; 
they have been calculated from the power spectra plotted 
in Fig. 2 by using (19). The first point to note is that ^e-i- 
differs from by less than 1% for angular scales larger 
than 3''; hence, the relative contribution caused by the 
source correlation is even smaller than that seen in the 
power spectra. This is due to the fact that the correlation 
function is a filtered version of the power spectra, how¬ 
ever with a very broad filter. This implies that even at 
small 0 scales, the correlation function is not dominated 
by large values of £, where the contribution from source 
clustering is largest, but low values of £ contribute signifi¬ 
cantly. The influence of source clustering on the ’ modes 
is larger, since the filtering function for those are narrower 
[i.e., J^x) is a more localized function that Jo(a::)], and ^e- 
differs from appreciably on scales below about 1'. 

Finally, in Fig. 4 we have plotted the aperture mea¬ 
sures. On scales below about 1', the dispersion of M± is 
larger than about 1% of that of Map. Hence, the ratio of 
these E- and B-mode aperture measures are very similar 
to that of the corresponding power spectra. 

The fact that Pce and Pb are very similar in amplitude 
means that by measuring Pb, one can make an approxi- 







Peter Schneider et al.: B-modes in cosmic shear from source redshift clustering 


11 



0[arcmin] 


Fig. 4. Aperture measures, for the same model as used in 
Fig. 2. Shown here is the dispersion of the aperture mass, 
(Mfp), the corresponding function in the absence of source 
correlations (noted by the subscript ‘0’) and which 

is the aperture measure for the B-mode. As expected from 
the power spectra shown in Fig. 2, and the fact that the 
aperture measures are a filtered version of the power spec¬ 
tra with a very narrow filter function, the B-mode aperture 
measure is considerably smaller than Map itself. 

mate correction of Pe, obtaining a value close to by 
subtracting Pe from Pe- Owing to the relative amplitude 
of these correlation-induced powers, such a correction may 
be needed in future high-precision measurements of the 
cosmic shear. 

5. Discussion 

We have shown that the clustering of galaxies from which 
the shear is measured leads to the presence of a B-mode 
in the cosmic shear field, in addition to providing an ad¬ 
ditional component to the E-mode. The reason for this 
effect in essence is the angular separation-dependent red- 
shift correlation of galaxies, which causes the mean of the 
product of the angular-diameter distance ratio along two 
lines-of-sight not to factorize, but to depend on 9. For a 
fiducial model considered in detail, the B-mode contribu¬ 
tion amounts to more than ~ 2% on angular scales below 
1' (or (. ^ 2.16 X 10^), and its relative importance quickly 
rises towards smaller angles. On substantially larger an¬ 
gular scales, however, the B-mode contribution is small. 
Furthermore, the additional E-mode contribution is very 
similar in size to the B-mode power, which will allow an 
approximate correction of the measured E-mode for this 
additional term. 


From an observational point-of-view, the most easily 
accessible quantities are the shear correlation functions 
^±, as one can easily deal with gaps in the data field. In 
Sect. 2 we have given explicit relations regarding how other 
two-point statistics of the shear can be calculated in terms 
of the shear correlation function. The finite support of 
the functions T± indicates that the aperture measures are 
more easily obtained from observational data than either 
the E- and B-mode correlation functions, or the E- and B- 
mode shear dispersions. Therefore, the aperture measures 
are the preferred method to check for the presence of a 
B-mode contribution in the shear data. 

We have varied some of the model parameters; in par¬ 
ticular, we have considered the case of lower mean source 
redshift (corresponding to a brighter flux threshold), and 
simultaneously increasing A(l'), such that the clustering 
length ro stays about the same. In this case we found 
a very similar ratio between the B- and E-mode power 
spectra as for the example considered in Sect. 4. We con¬ 
sider it unlikely that the observed B-mode in the present 
day data sets is due to the source clustering effect. The 
B-mode found in van Waerbeke et al. (2001) and Pen et 
al. (2002) can actually be used to search quantitatively 
for residual systematics. Its detection in van Waerbeke 
et al. (2001) was done by obtaining Map and M± by lay¬ 
ing down a grid of circular apertures on the data field. 
A more accurate measurement of (M^p) and (M^) has 
been obtained from the same data by Pen et al. (2002), 
by calculating them from the observed correlation func¬ 
tions , as in (HI) . In fact a subsequent analysis revealed 
that the B-mode measured in these data were essentially 
residual systematics caused by an overcorrection of the 
PSF, and can be corrected for (van Waerbeke et al. 2002). 
In this latter analysis, no significant B-modes are detected 
at small angular scales, but on scales above ^ 10', slightly 
significant values of M_l are detected; the effect discussed 
in this paper can certainly not account for them. 

The effect considered here seems to have been over¬ 
looked hitherto. Bernardeau (1998) considered the effects 
of source clustering on cosmic shear statistics and con¬ 
cluded that this source clustering can strongly affect the 
skewness and kurtosis of the cosmic shear, but to first or¬ 
der leaves the shear dispersion (and thus the power spec¬ 
trum) unaffected. Hamana et al. (2002) studied this ef¬ 
fect with ray tracing simulations, again concentrating on 
the skewness. Most of the other ray tracing simulations 
of weak lensing (e.g., van Waerbeke et al. 1999; Jain et 
al. 2000) assnmed all sources to be at the same redshift, 
in which case the additional power discussed here does 
not occur. Lombardi et al. (2002) calculated the effect of 
source clustering on the noise of weak lensing mass maps, 
showing that it can provide a significant noise contribution 
in the inner regions of clusters. 

It must be pointed out that the effect considered here 
is unrelated to other lensing effects which in principle 
could generate a B-mode, such as lens-lens coupling or the 
break-down of the Born approximation (see Bernardeau et 
al. 1997 and SvWJK for a discussion of these two effects on 
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the skewness). Numerical estimates (e.g., Jain et al. 2000) 
show that these latter two effects are very weak. Bertin 
& Lombardi (2001) considered the situation of lensing by 
two mass concentrations along the line-of-sight, where a B- 
mode is generated by a strong lens-lens coupling, but the 
fraction of lines-of-sight where this occurs is tiny. Another 
effect which could in principle generate a B-mode from 
lensing is the fact that the observable is not the shear it¬ 
self, but the reduced shear (Schneider & Seitz 1995). In 
the appendix we show that this effect is indeed negligible. 

Like the intrinsic alignment of galaxies, which can yield 
a spurious contribution to the measured cosmic shear, the 
source clustering effect can in principle be avoided if red- 
shift estimates of the source galaxies are available. In that 
case, by estimating the shear correlation function, pairs of 
galaxies with a large likelihood to be at the same distance 
can be neglected. In contrast to the intrinsic correlation of 
galaxies, the B-mode from source clustering appears to be 
fairly insensitive to the redshift distribution of the source 
galaxies, provided the clustering length is kept fixed. 


Appendix A: B-mode from reduced shear? 

The shear is not directly an observable, but is estimated 
from the image ellipticities of distant galaxies. The expec¬ 
tation value of the image ellipticity, however, is not the 
shear, but the reduced shear g = 7/(1 — k). Hence, the 
correlation of the observed ellipticities is the correlation 
of the reduced shear, not the shear itself. In cosmic shear, 
|k| <?; 1 nearly everywhere, and so the difference between 
shear and reduced shear shall not play a big role. However, 
at least a priori, this effect cannot be neglected, as seen 
from the following argument: 

The skewness S 3 = where X is a mea¬ 

sure of shear (such as Map, or the reconstructed k) has 
been calculated by van Waerbeke et al. (2001) to be of or¬ 
der a few hundred. On a scale of about one arcminute, 
~ 5 X 10“"'^, so that (A^) ~ 0.1 (A^), taking 
S 3 ~ 200 for the top-hat smoothed k. The difference be¬ 
tween the correlation functions involving g and those in¬ 
volving 7 is in principle of the same order-of-magnitude 
as and thus can be present at the level of a few per¬ 

cent, and there is no reason why it should not contain a 
B-mode contribution. 

We define the correlation functions 

C±( 6 ') = {gtiO)gt{e)) ± {gx{0)gx{0)) (A.l) 

and choose 9 = {9, 0), so that gt = —gi, gx = — 92 - Using 
the approximation g « 7(1 -|- k), valid for |k| <C 1 , we 
obtain 


Replacing the shear and convergence by their Fourier 
transforms, this becomes 

H/H/H 

X -b cos [2 (/3i t 02 )] , (A.4) 

where, as before, 0i is the polar angle of £i. The triple cor¬ 
relator vanishes unless the sum of the wave-vectors equals 
zero; one defines the bispectrum by 


{k{£i)k{e2)k{i3)) = {27 t)^ S{£i+£2+£3)b{£i,£2,£3) .(A.5) 
Performing the .^ 3 -integration in ( 0 ) yields 

X -b cos [2 (/3i =f 02)] ■ (A. 6 ) 


The function b{£i,£ 2 , —£i — £ 2 ) has three independent ar¬ 
guments, namely the moduli ii and £ 2 , and the angle 
0 = 01 — 02 between the two £-vectors. We therefore write 
b{£i,£ 2 , —£i —£ 2 )= &(£i, £ 2 , 0), make use of the symmetry 
in the integrand of o, and replace the /3i-integration 
by one over 0 : 


Ae+(d) = 2 


^0 

n27r 


dil h 


d£. 


2 ^2 


/O 


(2^)2 


d0 b{£i,£2,0) 


>0 


>0 


d 02 cos( 2 (/)) 

1 r 

TT 


/O 

/*27r 


X / d(/) cos(2(()) 6 (£i,£ 2 , 0) ; 

Jo 

analogously, one obtains 


(A.7) 


A^e) = - 

TT 


d£i £1 

( 2 ^) 


J4(£i0) 


d £2 


( 2 ^) 


d0 cos{ 20 )b{ii,£ 2 , 0 ) ■ 


(A. 8 ) 


Inserting these expressions into ( |l 8 | ) and making use of 
( 0 ) , one immediately sees that the reduced shear does not 
yield any B-mode contribution, and that the correlation 
functions for the reduced shear are 


e±io) = [ [p.{£)+p^^Hi) 


(A.9) 


where 


^^(0) = e±(0) + A^±(0), 

where 

X^±{9) = ( 71 ( 0 ) 71 ( 0 ) [k( 0 )-b k( 0 )]) 

± ( 72 ( 0 ) 72 ( 0 ) [k( 0 )-b k( 0 )]) . 


(A.^) 


fdH' 

■ 2 (£-£')^ 

j (271)2 

[|£| 2|£'|2 J 


&(£,£',-£-£').(A.10) 


Thus, considering the reduced shear yields an additional 
E-mode power to the one obtained from considering the 
(A.3) shear itself. 
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